多排叶片内部流动以及叶片排与排气系统相互作用研究

本文采用高分辨率格式和多块多网格方法求解雷诺平均Navier-Stokes方程,数值模拟多排叶片内的三维粘性流动;使用数值激盘模拟叶片排的存在,研究叶片排与非轴对称排气部件之间的相互作用以及复杂的内部流动。文中描述了数值方法,给出了NASA透平导叶和单级透平内部流场的数值结果及其与相应实验结果的对比,也给出了多级透平内部流场的数值结果,以及透平与排气部件之间的耦合流场的结果。     

热动力学问题的数值计算

 对含热传导的流体动力学方程组，用有限元方法进行数值求解。采用傅里叶热传导计算热流、用热流连续条件计算单元间接触面的温度、用三角形传热法计算体单元表面的热流，考虑各向同性弹塑性流体材料模型、三项式物态方程和导热系数与状态的相关性，给出了傅里叶热传导、接触传热、热应力应变效应、以及混合物冲击压缩特性等算例。对混合物冲击温度的数值模拟表明，小颗粒混合物在冲击压缩过程中，颗粒间的温度有差别、稍有波动，并随时间趋向于一致，以至热平衡。     

The influence of temperature gradients on the kinetic boundary layer problem for a condensing droplet

We extend an earlier method for solving kinetic boundary layer problems to the case of particles moving in aspatially inhomogeneous background. The method is developed for a gas mixture containing a supersaturated vapor and a light carrier gas from which a small droplet condenses. The release of heat of condensation causes a temperature difference between droplet and gas in the quasistationary state; the kinetic equation describing the vapor is the stationary Klein-Kramers equation for Brownian particles diffusing in a temperature gradient. By means of an expansion in Burnett functions, this equation is transformed into a set of coupled algebrodifferential equations. By numerical integration we construct fundamental solutions of this equation that are subsequently combined linearly to fulfill appropriate mesoscopic boundary conditions for particles leaving the droplet surface. In view of the intrinsic numerical instability of the system of equations, a novel procedure is developed to remove the admixture of fast growing solutions to the solutions of interest. The procedure is tested for a few model problems and then applied to a slightly simplified condensation problem with parameters corresponding to the condensation of mercury in a background of neon. The effects of thermal gradients and thermodiffusion on the growth rate of the droplet are small (of the order of 1%), but well outside of the margin of error of the method.     

自由电子激光振荡器的二维数值模拟

在柱二维坐标系中,用数值计算方法求解了摇摆器内电子模拟方程组、光场方程组和无源区(不含摇摆器区)经非自适应方法变换后的无源旁轴波动方程并设计了相应的程序R_2D。用该程序计算的数值结果表明,在有源条件下二维数值模拟结果与高斯光束近似下的数值模拟结果符合较好;取完全相同的计算条件R-2D程序与洛斯阿拉莫斯实验室(LANL)的FELEX程序的计算结果符合较好;初步认为LANL实验数据的物理图象合理,数值结果基本可信;最后给出北京自由电子激光器(BFEL)的一些物理数据。     

A parallel inexact newton method for stochastic programs with recourse

A parallel inexact Newton method with a line search is proposed for two-stage quadratic stochastic programs with recourse. A lattice rule is used for the numerical evaluation of multi-dimensional integrals, and a parallel iterative method is used to solve the quadratic programming subproblems. Although the objective only has a locally Lipschitz gradient, global convergence and local superlinear convergence of the method are established. Furthermore, the method provides an error estimate which does not require much extra computation. The performance of the method is illustrated on a CM5 parallel computer.This work was supported by the Australian Research Council and the numerical experiments were done on the Sydney Regional Centre for Parallel Computing CM5.     

曙光一号自由电子激光器的理论计算

系统总结曙光一号自由电子激光器理论计算的主要结果:包括曙光一号装置主要参数的选取和理解;磁场失谐曲线的计算;常参数摇摆器和变参数摇摆器的主要结果;高阶波导模的贡献;电子束参数扰动对激光性能的影响;空间电荷效应等。计算结果表明,常参数摇摆器激光输出功率可达80MW,效率约50％;变参数摇摆器激光输出功率可达250MW左右,效率约16％。     

Numerical solution of eigenvalue problems for Hamiltonian systems

This paper discusses the numerical solution of eigenvalue problems for Hamiltonian systems of ordinary differential equations. Two new codes are presented which incorporate the algorithms described here; to the best of the author’s knowledge, these are the first codes capable of solving numerically such general eigenvalue problems. One of these implements a new new method of solving a differential equation whose solution is a unitary matrix. Both codes are fully documented and are written inPfort-verifiedFortran 77, and will be available in netlib/aicm/sl11f and netlib/aicm/sl12f.     

Numerical propagation of partially coherent laser beams through optical systems

The propagation of partially coherent beams through optical systems is computed numerically in one transverse dimension. The optical system is divided into different elementary segments, through which the propagation of light can be calculated by appropriate operators, working on the coherence function or the Wigner distribution function respectively. For the necessary changes between these two functions describing the partially coherent beams, the use of the remarkable z-transform seems to be an advantage. With this algorithm the grid and the resolution in the spatial frequency domain can be arbitrarily chosen in contrast to the usual Fourier transform, the influence of phase aberrations on the focusability of Gauss-Schell model beams is discussed as an application example of the numerical model. With the help of this tool, practical beam guiding systems can be simulated for use with multimode laser radiation.     

斜孔气膜冷却数值模拟分析

本文通过数值模拟分析圆柱孔和扩散孔的单斜孔气膜冷却特性,考察复合角、孔型和吹风比对流场和气膜冷却效果的影响。结果表明,复合角的引入使气膜侧向分布更宽,但冷却效果沿主流方向衰减更快。适中吹风比得到的气膜能更有效的保护壁面。在相同吹风比和复合角条件下,扩散孔的气膜冷却效率比圆柱孔更好,且冷却更为均匀持久。     

Numerical Computation of Figure-eight Solution for 3-body Problems

The main goal of this paper is to compute the Figure-eight solutions for the planar Newtonian 3-body problem with equal masses by finding the critical points of the functional associated with the motion equations of 3-body in plane R2. The algorithm adopted here is the steepest descent method, which is simple but very valid for our problem.